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Abstract 



A statistical physics model for the time evolutions of stock portfolios is 
proposed. In this model the time series of price changes are coded into the se- 
quences of up and down spins. The Hamiltonian of the system is introduced 
and is expressed by spin-spin interactions as in spin glass models of disor- 
dered magnetic systems. The interaction coemcients between two stocks are 
determined by empirical data coded into up and down spin sequences using 
huctuation-response theorem. Monte Carlo simulations are performed and 
the resultant probability densities of the system energy and magnetization 
show good agreement with empirical data. 
The data analysis and modeling of nnancial markets have been hot research subjects for 
physicists as well as economists and mathematicians in recent years p|. The non-Gaussian 
property of the probability distributions of price changes, in stock markets and foreign 
exchange markets, has been one of main problems in this field [0] [^J. From the analysis 
of the high-frequency time series of market indices, e.g., S&P500, Nikkei225, a universal 
property was found in the probability distributions. The central part of the distribution 
agrees well with Levy stable distribution 01], while the tail deviate from it and shows another 
power law asymptotic behavior. The scaling property on the sampling time interval of data 
is also well described by the crossover of the two distributions. Several stochastic models of 
the nuctuation dynamics of stock prices are proposed, which reproduce power law behavior 
of the probability density 0. The auto-correlation of nnancial time series is also an 
important problem for markets. There is no time correlation of price changes in daily scale, 
while from more detailed data analysis an exponential decay with a characteristic time r = 4 
minutes was found 0. The fact that there is no auto-correlation in daily scale is not 
equal to the independence of the time series in the scale. In fact there is auto-correlation 
of volatility (absolute value of price change) with a power law tail ||. Recently, the 
cross-correlation between pairs of stock issues was deeply investigated for the time series 
of price changes, and the hierarchical structure in the subdominant ultrametric space 
was found ||. This result suggests the complex collective time evolution of nnancial 
markets as in frustrated disordered systems like spin glass || in which the ultrametricity 
has been established. Those problems listed here have been studied using the methods and 
the concepts developed in material sciences especially in the studies of complex systems. 

In this paper, a statistical physics model for the collective time evolutions of stock port- 
folios is proposed. Portfolio is a set of stock issues. In this model we deal with the time 
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series of price changes coded into the sequences of up and down spins. A sample of coding 
procedure is shown in Fig. |l[ The Hamiltonian of the system is introduced and is expressed 
by spin-spin interactions as in spin glass models of disordered magnetic systems. The inter- 
action coemcients between two stocks are phenomenologically determined by empirical data. 
They are derived from the covariance of sequences of up and down spins using Auctuation- 
response theorem. We investigate the stocks listed in Dow-Jones industrial average as a 
portfolio for the test of this model. Monte Carlo simulations using Gibbs weight as a tran- 
sition probability reproduce the probability densities of the energy and the magnetization 
of the portfolio, whose dehnitions are given later. 

We start with the Hamiltonian expression of our system that contain N stock issues. It 
is a function of the conhguration S consisting of N coded price changes S{ (i = 1, 2, N) 
at equal trading time. The interaction coemcients are also dynamical variables, because the 
interactions between stocks are thought to change from time to time. We divide a coemcient 
into two parts, the constant part J^, which will be phenomenologically determined later, 
and the dynamical part bJ^. The Hamiltonian including the interaction with external helds 
hi (i = 1, 2, N) is denned as 

SJ 2 - 

H[S, SJ, h]=J2[—^- (J Z3 + SJ^SiSj)] - £ hiSi. (1) 

<i,j> ZL * l 3 i 

The summation is taken over all pairs of stock issues. This form of Hamiltonian is that 
of annealed spin glass |§. The Auctuations SJij are assumed to distribute according to 
Gaussian function. A statement referring to the justihcation of this assumption will be 
made later (see Fig. ^). The main part of statistical physics is the evaluation of partition 
function that is given by the following functional in this case 

Z[h) = £ / II ( Ja\/2 eM-H[S,SJ,h\). (2) 
{Si} <i,j> y^^o) 

The integration over the variables SJij is easily performed and gives 

Z[h]=Aj2eM-H ef f[S,h}). (3) 

{Si} 

Here the effective Hamiltonian H e ff[S, h] is dehned as 

H ef f[S,h]=— £ JijSiSj — ^hiSi, (4) 

<i,j> i 

and A = exp(l/2 Ajj) is just a normalization factor which is irrelevant to the following 
step. This form of Hamiltonian with constant Jij is that of quenched spin glass 0. 

The constant interaction coemcients Jij are still undetermined. We use nuctuation- 
response theorem which relates the susceptibility \ij with the covariance C^- between dy- 
namical variables in order to determine those constants, which is given by the equation 

= dj. (5) 

h=0 



drrii 
Xij = dh~ 
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The notations in this paper are m^ =< Si > and Cy =< Si — m^ >< Sj — rrij >. For 
the evaluation of the left side of the equation (^), we use Thouless-Anderson-Palmer (TAP) 
equation for quenched spin glass |10| 



rrii 



tanh(^ Jijmj + J?j(l ~ rn))rrii)- (6) 



We will neglect the third term of self-reaction |TTJ and the nonlinear terms ~ 0(mf) |T2| in 



the equations @. The equation (||) and the linear approximation of the equation (|(J) yield 
the equation 

y"X^ifc — JikjCkj = 5ij. (7) 
k 

Interpreting C^ as the time average of empirical data over a observation time rather than 
ensemble average, the constant interaction coemcients is phenomenologically determined 
by the equation (^). 

We investigate two data sets of a portfolio containing iV = 30 stock issues listed in Dow- 
Jones industrial average for the test of this model. One is the time series in the period from 
16-May-2000 to 21-Jun-2000 and the other is in the period from 15-Aug-2000 to 26-Oct- 
2000. The number of data amounts to about 30, 000 x 30 = 900, 000. The sampling interval 
is 1-minute, i.e., the interval two successive time stamps of samples are 1-minute. The time 
series of price changes are coded into the sequences of up and down spins (see Fig. |TJ again). 
This coding is thought to be a kind of coarsening. The covariance C^ for 435 pairs of i and 
j are derived from the coded data. Then the interaction coemcients J^ are calculated by 
the equation (|7|). The energy spectra of the system, simply the portfolio energy, is denned 
as the eigenvalues of the Hamiltonian H e f f [S, 0] . The probability density of the portfolio 
energy can be obtained in two ways. We can calculate the probability density from data by 
the equation 

ApJ AE 
p(E)AE = P( E - — <E<E + —). (8) 

The results for two data sets are shown in Fig. [| (O) with the results of Monte Carlo 
simulations (□) that will be explained below. The use of the theoretical equation 

p(E)AE = n(E)^ (9) 

is the other way. The function n(E) is the density of states, which can be numerically 
obtained by random sampling from 2 30 conhgurations. The comparison between the two 
results is given in Fig. |3|. The theoretical lines explain well the empirical data except the 
ranges E < —4 and 1 < E of small numbers of events, indicating the canonical distribution 
of the system conngurations S as P(S) = exp(—H e ff[S,0])/Z[0]. 

For another test of our model, Monte Carlo simulations using Gibbs weight as a transition 
probability are performed. The transition probability from a connguration S to 5" is given 
as 

exp(-5H) if 5H>0 
W (S -+S') = { (10) 
1 otherwise, 
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where 5H = H e ff[S', 0] — H e ff[S, 0]. Monte Carlo simulations for the probability density of 
the portfolio energy and the system magnetization m = (1/N) Yli=i S{ are given in Fig. |2| 
(□) and Fig. [| (□) respectively. They also show good agreement with empirical data. For 
the final consistency check of our model, we calculate the covariance C™ between pairs of 
two stocks by Monte Carlo simulation. The frequency distribution of 5 = {(Cy — C™ c )/Cij\ 
is given in Fig. |5], which indicate the reproduction of C^ to the extent of our approximation. 

Here I want to refer to the justincation of the first term in the equation ([!]), which 
indicates that the Auctuations 5Jij distribute according to Gaussian function. The direct 
verification of this assumption is impossible, because 5Jij change from time to time. In- 
stead if we take the observation time long enough, the justincation of this assumption is 
mathematically obtained by central limit theorem. For the purpose of the test whether our 
observation time is long enough, we put the following test to the data in the period from 
15-Aug-2000 to 26-Oct-2000. First we take successive 50 trading days from the data, and 
then divide it into 10 pieces. Next we derive j£ ( k = 1, 10) for each piece of data to 
obtain the distribution of 435 x 10 5J^ = J^ — Jij. The result is shown in Fig. || supporting 
the statement of central limit theorem. 

In this paper, we gave a fully consistent phenomenological model for stock portfolios, 
which is expressed by the effective Hamiltonian (|j). This model will be also applicable to 
other nnancial markets that show collective time evolutions, e.g., foreign exchange market, 
options markets, inter-market interactions. In our model, however, J^ must be determined 
by data via the calculation of Cy . The problem how we can go beyond this phenomenological 
standing point is unsolved. There may be a hint in the universal property of the covariance 
matrbc C = [Cy] [|14j or (and) the existence of the correlation between a pair of N(N — l)/2 
Jij corresponding to two separated time periods investigated above |13 . 
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FIGURES 



FIG. 1. A sample of the coding procedures for a virtual data. The nrst column is the serial 
number of data sampling, in which means the time origin of observation. Samplings go along 
trading time at a fixed interval, e.g., sampling every 1-minute. From the price of i-th sample Y(i), 
the i-th price change is dehned as Z(i) = Y(i) — Y(i — 1). The i-th code S(i) = 1 (if Z(i) > 0), 
S(i) = -1 (if Z(i) < 0) and S(i) = S(i - 1) (if Z(i) = 0). 

FIG. 2. The probability density p(E) of portfolio energy E. : Fhe empirical probability 
density (AE = 0.1 in the equation (§)). Monte Calro simulation of 2 17 samples, which is 
explained later. (a) The result for 9994 x 30 coded data of the time series sampled at 1-minute 
time inteiwal in the period from 16-May-2000 to 21-Jun-2000. (b) The result for 19992 x 30 
coded data of the time series sampled at 1-minute time interval in the period from 15-Aug-2000 to 
26-Oct-2000. 

FIG. 3. Semi-log plot of the probability weight p(E)n(0)/p(0)n(E) of portfolio energy E. The 
probability density p(E) is the empirical result shown in Fig. 2, and the density of states n(E) is 
numerically obtained by 2 17 random sampling from 2 30 conhgurations. : The empirical proba- 
bility weight. Solid line: Gibbs weight e~ E (no ntting parameter). (a) The result for the same 
data as in Fig. 2(a). (b) The result for the same data as in Fig. 2(b). 

FIG. 4. The probability density of the system magnetization dehned in text for the coded 
data in the period from 15-Aug-2000 to 26-Oct-2000. O : The probability density whose mean 
\x = —0.0134, variance a 2 = 0.0744. Dashed line: Gaussian hmction with the same mean and 
variance as of the probability density. □: Monte Calro simulation of 2 19 samples 

FIG. 5. The frequency distribution of 5 dehned in text. 

FIG. 6. The probability density of the distribution of 435 x 10 5 whose dennition is given in 
text. : The probability density whose mean /i = —0.000324, variance a 2 = 0.000679, skewness 
fj\ = —0.0467 and kurtosis /?2 = 0.0604. Dashed line: Gaussian hmction with the same mean and 
variance as of the probability density. 
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price change 
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81.7500 
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0.0625 
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0.1250 
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4 


81.8125 
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81.8125 
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-1 


6 


81.7500 


-0.0625 


-1 
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< S < 0.1 


227 


52.18 


0.1 < S < 0.2 


113 


78.16 


0.2 < S < 0.3 


56 


91.03 


0.3 < S < 0.4 


19 


95.40 


0.4 < S < 0.5 


8 


97.24 


0.5 < S < 0.6 


7 


98.85 


0.6 < S < 0.7 


2 


99.31 


0.7 < S < 0.8 


2 


99.77 


0.8 < S 


1 


100 



